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1 1 Introduction 

2 A central goal of flavor physics is to measure the angle 7 = arg (~ y^~ ) in the 

3 Cabibbo-Kobayashi-Maskawa (CKM) [U [2] mixing matrix, which is currently known 

4 to a precision of about 10-12° |3J. The theoretically cleanest methods employ B — > 

5 DK decays, where the sensitivity to 7 results from the interference between b — > c and 
e b — y u transitions. Since both transitions are 0(X 3 ) in the Wolfenstein parameter [1], 
7 large CP violating asymmetries are expected. One powerful class of methods utilize 
s B~ — > DK~ where the D is detected in either a CP eigenstate [Zj, a flavor-specific 

9 mode |6], or a multi-body decay [!8J. An advantage of these decays is that they do 

10 not require knowledge of the 6-hadron flavor at production (flavor tagging), and only 

11 rely on measuring the time integrated rates. Another powerful method to extract 

12 7 is to perform a time-dependent analysis of B Q S — > DfK~ [9j [101 E] and B® — > 

13 D^K^n + n^ . Time-dependent analyses of B® — > D+K~(ty + 'k~) are only possible at 
u hadron colliders, and are a unique capability of LHCb. 

15 The time-dependent decay rates of B° s and B° s to a flavor-specific final state, 
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where A/ is the decay amplitude A(B® -> /) and A/ = (q/p)(A f /A f ) = |A/|e i(A - (7 ~ 2/3s)) . 
Here, |A/| and A are the relative magnitude and strong phase difference between the 
b — y u and b — > c transitions, and 2/3 s is the phase of mixing. The complex coeffi- 
cients p and q relate the B° s meson mass eigenstates, Bh,l, to the flavor eigenstates, 
B® and Tf s via: 

Bl = pB° s + qB° s 
B H = pB°-qB* ' 



IpI 2 + M 2 



i . 



(3) 



22 Similar equations can be written for the CP-conjugate decays, replacing Af by Aj = 

23 A{W S -> /), A/ by A 7 = (p/q)(Aj/Aj), \p/q\* by \q/p\\ C f by Cj, S f by Sj, and D f 

24 by Dj. The CP asymmetry observables Cf, Sf, Df, Cj, Sj and Dj are then given by 
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25 Since CP violation in mixing is expected to be below the percent level, it follows that 

26 \q/p\ = 1, |A/| = |Aj|, and consequently Cf = Cj. Thus there are five observables 

27 that depend on the 3 physics parameters of interest: |A/|, A and 7 — 2(3 S . Similar 

28 expressions are applicable to B® — >■ DfK^n + TT~, however, there is a potential dilution 

29 due to the varying strong phase across the D+K~n + 'K~ Dalitz plane. 

30 In this article, we present the first measurements of these five CP observables. 

31 First observations of the B° s -> D+K-ir+n-, B° ->■ D+R-^n- and B s ->■ D sl (2536) + 

32 decays are also presented, along with measurements of their relative branching frac- 

33 tions. All results are based on 1.0 fb _1 of integrated luminosity recorded in 2011 

34 by the LHCb experiment. More detailed documentation of the B® — > DfK~ and 
35 
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analyses can be found in Refs. p2] and [15] . respectively. 
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36 2 Event Selection 



37 Signal candidates are formed by reconstructing Df — > K + K 7r + , Df — > n + n n + 

3B and Df K + n~n + . For the B ^ -)■ D+K^n + n^ and 5° ->■ Dfn^n + 7i^ candidates, 

39 only the Df — > K + K~n + decay is considered. The Df candidates are required to 

40 form a good quality vertex, be spacially well separated from any primary vertex (PV), 

41 and have an invariant mass consistent with the known mass (within about 3 times 

42 the mass resolution). Multivariate selection algorithms are employed to suppress the 

43 combinatorial background, and typically have a signal efficiency of 80-90% while re- 

44 jecting about 85% of the combinatorial background. Invariant mass distributions for 

45 candidates are shown in Fig.[TJfor the higher signal yield B° s — > D+ir~ decay, show- 
ing that clean signals are achievable even in the suppressed -D+ decay modes. Tighter 
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Figure 1: Invariant mass distributions for Df candidates in the B° s — > D+7r~ data 
sample, for (left) K + K~tt + , (middle) K + n~ir + , and (right) tt + 7t~tt + final states. 

46 

47 particle identification requirements are applied to the K~ or K~ix + -k~ recoiling from 

48 the to suppress cross-feed from the favored B° s — > D+tt~ and B° s — > Df-n-~7T + 7T~ 

49 decays. For the B® — > D+7t~7t + tt~ and B° s — > D^K^-k + -k~ decays, the invariant mass 

50 of the 7r~7r + /T~ and K^-k + -k^ systems are restricted to be below 3000 MeV/c 2 . 

n 3 Analysis of B° s -> D+tT and B° s -> D+K~ 

52 The invariant mass distributions for B® — > Dfn~ and B° s — > D+K~ are shown 

53 in Figs. [2] and [3j All three decay modes have approximately equal B° s mass 

54 resolutions, and are summed together in these distributions. The signal shape is 

55 modeled as the sum of two Crystal Ball [2] functions, with one exponential tail on 

56 each side of the B° s signal peak. A number of specific backgrounds, due to either a 

57 missed particle (e.g. B° s — > D+p~, with the ir° undetected), a misidentified particle 

58 (e.g. B° s ->■ D+n~ reconstructed as B® D+K~), or both (e.g. B° s -D+p~ 
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59 reconstructed as B° s — > D+K ) are accounted for using either data or simulation 

60 to model the shape of these backgrounds. From an unbinned extended maximum 

61 likelihood fit, 27, 965 ± 395 B° s -> D+tT and 1390 ± 98 B° s ->■ D+K~ signal events 
are selected. The CP parameters are obtained by a fit to the decay time distribution 
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Figure 2: Invariant mass distributions B° s — > D+n~ candidates. The signal component 
is indicated by the dashed curve, and the backgrounds are indicated by the various 
color- filled (shaded, in B/W) curves. 

62 

63 of the B° s — > D^K~ signal candidates. Two methods have been developed. The first, 

64 referred to as sFit, uses sWeights [T5] obtained from the B Q S — > D+K~ mass fit to 
es statistically subtract the background contribution. The second method, referred to 
ee as cFit, is a conventional two-dimensional fit to the reconstructed mass and decay 
67 time. The advantage of the first method is that there is no need to model the time 
ee distribution of all the backgrounds, as they are statistically removed via the s Weights. 

69 The statistical subtraction, as presented here, uses events in the full mass fit region, 

70 and the subtraction of this background leads to a larger statistical uncertainty than 

71 if just a narrow signal region is used. For this reason, the second method is expected 

72 to give a smaller statistical uncertainty; however it requires an accurate model of 

73 the time distributions of the backgrounds that enter into the signal region. For the 

74 analysis presented here, the sFit provides the nominal result, and the cFit is used as 

75 a cross-check. 

76 The measurement of the CP parameters in B° s — > D^K~ requires a fit to the time- 

77 dependent decay rates. The fit accounts for (i) the acceptance versus reconstructed 

78 decay time, (ii) the decay time resolution, and (iii) the effective tagging efficiency The 
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Figure 3: Invariant mass distributions B° s — y DfK~ candidates. The signal com- 
ponent is indicated by the dashed curve, and the backgrounds are indicated by the 
various color-filled (shaded, in B/W) curves. 



79 functional form of the acceptance function is determined from simulated B° s — y D+tt~ , 

so and its parameters are determined in a fit to B® — y D^ir~ data, where the B° s lifetime 

si and mixing frequency, Am s , are fixed to 1.51 ps and 17.69 ps _1 [17], respectively. 

82 The average decay time resolution is about 50 fs, and is modeled by the sum of 

83 three Gaussian functions, whose parameters are determined from simulation. The 

84 Gaussian width parameters obtained from simulation are scaled up by 1.15 to account 
as for better resolution in the simulation than in data; this factor is determined by 
se comparing the width of the zero decay time component of prompt plus one random 
87 track in data and simulation. For the flavor tagging, only opposite side (OS) taggers 
as are currently used. These algorithms exploit the correlation in flavor between the 

89 signal b hadron at production, and the other b hadron in the event (referred to as 

90 the tag-fe). In particular, the charge of either an electron, a muon, or a kaon that 

91 does not come from any pp interaction vertex (or the signal b), or the charge of 

92 another secondary vertex in the event, provide information on the flavor of the tag-6 

93 hadron. Because bb are produced in pairs, this translates into a flavor determination 

94 of the signal B®. The OS flavor tagging algorithm was initially tuned using simulated 

95 decays, and then re-optimized and calibrated to obtain the largest effective tagging 

96 efficiency using the self-tagging B + — y J/ipK + and B° —y D*~fi + v decays in data. In 

97 general, the performance of the OS tagging algorithms are independent of the signal 

98 6-hadron decay, and have a combined effective tagging efficiency of eD 2 = 1.90% for 
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99 B° s — > D+K . Further details of the tagging algorithms can be found in Ref [15] . 

100 In the fit to B° s — > DfK~, the following parameters are fixed: Am s = 17.69 ps _1 , 

101 t_ Bs = 1-51 ps and AT S = T SiL - T S , H = 0.105 ps" 1 [TF]. About 60% of the 

102 B® — > DfK~ candidates have no flavor tag; the time-dependent decay rates for 

103 these untagged decays is given by the sum of the two expressions in Eq. El and 

104 the sensitivity to 7 enters through the hyperbolic sine term. The decay time dis- 

105 tribution of B® — > DfK~ signal decays and projections of the fitted are shown in 

106 Fig. HI The projections show the four possible tagged decays, B® — > DfK T and 

107 B® — > DfK^, as well as the untagged time-dependent decay rates (B®, 5°) — > D~K + 
and (B®,B®) — > D+K~. The fitted values for the CP parameters are 




Figure 4: Distribution of reconstruct decay time for B° s — > DfK~ signal decays 
(points with error bars), along with the results of the fit. Projections of the decay 
rates versus the decay time for the four possible flavor tagged decays, and the two 
untagged decays. 
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108 



1.01 ±0.50 ±0.23, 
1.25 ±0.56 ±0.24, 
0.08 ±0.68 ±0.28, 
1.33 ±0.60 ±0.26, 
0.81 ±0.56 ±0.26, 

loo where the first uncertainties are statistical and the second are systematic. Several 
no sources of systematic uncertainty have been considered. The dominant sources are 
in due to the precision on the effective flavor tagging efficiency (0.16cr st at-0.23cr s tat), vari- 
n2 ations in the parameters that are fixed in the default fits (0.15<7 s tat-0.42cr st at), and 
n3 the correlation between the mass of specific backgrounds and their reconstructed de- 
H4 cay time (0.08cr sta t-0.27cr s tat), where these uncertainties are expressed as a fraction 
us of the statistical error. These are the first measurements of the CP parameters in 
lie B° s — > DfK~. With additional data and analysis refinements, reduction in both the 
ii7 statistical and systematic uncertainties are expected. 

us 4 First observation of B s — > DfK~7T + 7T~ and B° s — > 
D s i(2536) + tt- 

120 The decay B° s — >• DfK~7T + 7T~ can be analyzed in a similar way to B° s — > D+K~ to 

121 measure the weak phase 7. While this decay has not yet been observed, if one uses 

122 B° and B~ decays as a guide, it would naively be expected that its branching fraction 

123 is 1.5-2.0 times larger than B® — > DfK~, making this a potentially attractive decay 

124 mode to explore. The first step in such an analysis is to firmly establish an observation 

125 of this decay and measure its branching fraction (here, relative to B° s — > D+n~n + n~). 

126 While searching for this decay, the decay B° — > D+K^n+n^ is also observed and its 

127 branching fraction is measured relative to B° s — > DfK~7T + 7T~ . 

128 With the previously defined selections, Fig. [5] shows the invariant mass distri- 

129 butions for (left) B° s — > Dfn~n + n~ candidates and (right) B°, ' s — > K~n + n~ 

130 candidates. Significant B° s signals are seen in both spectra, and a B° signal is seen 

131 in the D^K~tt + tt~ mass distribution. The main sources of background are B° s — > 

132 D*+it-it + it- (to B° s ->- D+n-iT+n-), and B° s ->- D+tt-tt+tt-, B° s ->- D* s + ti + , 

133 and B ^ — > DI + K~tt + tt~ (to S9s — > D+K~7t + tt~). Their shapes are taken from 

134 simulation, with parameters that are allowed to vary within their uncertainties. 

135 Yields of 5683 ± 83 B° s D+tt-tt+tt", 216 ± 21 B° s D+K~-K + 7r- and 402 ± 33 
we B° — > K^-k + -k^ are observed. After correcting for the relative efficiencies, the ratio 
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137 of branching fractions are measured to be 

= (5.2 ±0.5 ±0.3) x 10 
= 0.54 ±0.07 ±0.07, 



B(B° S -4 D+K-ti+ti-, 

£(5° -4 D+TT-TT+TT-; 
£(5° -4 D+K-TT+iT- 



B(B° S -4 D+K-TT+ 
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where the uncertainties are statistical and systematic, respectively. These are the 
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Figure 5: Invariant mass distribution for (left) B Q S -4 D+vt~/t + /t~ candidates and 
(right) B9 x -4 DfK~ir + ir~ candidates. The fitted signal (dashed lines) and back- 
ground shapes (shaded/hatched regions) are shown, as described in the text. 



138 

139 first observations of these decays. Since B° s -4 D+tt^tt+tt^ has a branching fraction 
wo that is about twice as large as B° s -4 Dfn~, and B(B° S -4 D+K~) ~ 0.09 x -4 

141 D+-R-) [18J, it follows that B(B° S -4 -D+?T~7r + /T~) is at least as large as B(B° S -4 

142 Dfn~), or as much as 50% larger. The i3(_B° -4 DfK~7T + 7T~) is also sizeable, and is 

143 likely dominated by contributions where an extra ss pair is produced in addition to 

144 the weak decay (see Ref. [33] for more details). 

145 The B° s -4 D+7t~7t + 7t~ decay has also been analyzed to search for intermediate 
us excited D S j states. For B° s -4 D+7t~7t + tt~ candidates within 40 MeV/c 2 of the B° s 
147 signal peak, the mass difference, AM = M(Dfn^7i + ) — M(D+) is computed for both 
us 7r~7r + mass combinations. The resulting mass difference spectrum is shown in Fig. [6j 
149 The signal is fit with a Breit-Wigner convolved with a Gaussian resolution function 
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whose width is fixed to the expected AM resolution. A signal of 20.0 ±5.1 events is 
observed with a AM value and width consistent with the D sl (253Q) + state. Applying 
corrections for the relative efficiency, the ratio of branching fractions is measured to 
be 

8(50^(2536)^-, Dt^Dj^) = ±lQ± x 1Q _ 3 
B(B° S ->• D+n-n+n-) 



154 The excess of events is 5.9 standard deviations over the expected background, thus 
establishing the first observation of this decay. 
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Figure 6: Distribution of the difference in invariant mass, M(D+n~n + ) — M(Df), 
using B° s — > D+7t~7t + 7t~ candidates within 40 MeV/c 2 of the known B° s mass (points) 
and in the upper B° s mass sidebands (filled histogram). The fit to the distribution is 
shown, as described in the text. 



156 5 Summary 

is? First measurements of the CP observables in the B° s — > DfK~ decay have been re- 

158 ported. With the larger data sample recorded in 2012, and the larger data set antici- 

159 pated in the future, this decay will contribute significantly to the determination of the 
wo weak phase 7. First observations of the B Q S — > K~n + n~ and B° — > K~n + n~ 
i6i are also reported. The former can be used in a similar way to B° s — > D+K~ to 
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162 extract 7. After including D+ — > and D+ — > K^n + n^ decays, and re- 

163 optimizing the selection for B® — > DfK~n + 'K~ only, the yield in this mode more 

164 than doubles with a comparable signal-to-background. The yield in this mode is 

165 therefore expected to have about 35-40% of that obtained in B° s — > DfK~. The 
lee B® — >■ /} s i(2536) + 7r~ decay is also observed for the first time, and its branching frac- 

167 tion relative to B° s — > _D+?T~?T + 7r~ is presented. 
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